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I. INTRODUCTION 


Recently, quantum information theoretic measures like 
decoherence [IHll, entanglement and fidelity have 
become the subject of immense interest. Close to the 
quantum critical point (QCP) of a quantum many-body 
system, these quantities show peculiar behaviors and 
hence can be regarded as a tool to detect the QCP 0. 
Decoherence that signifies the loss of coherence of the 
system when it interacts with the environment is an im¬ 
portant observable for quantum computation. 

In this context, Loschmidt echo (LE), which quanti- 
hes the decoherence is studied extensively HQ. The 
LE is defined as the square of overlap of the two wave 
functions |'0(i)) and |'0o(i)) evolving with two different 
Hamiltonians H and Hq, respectively, z.e., 

= ( 1 ) 

Initially, both the states are prepared in the ground state 
of Hq. LE provides information about how small pertur¬ 
bations during an evolution can result to the decoher¬ 
ence of the state of the system, thus being an important 
quantity for information processing and storage. On the 
other hand, LE can also be used to detect the presence 
of a QCP by showing a sharp dip at the QCP of the 
Hamiltonian when H and Hq are close to each other. 

At the same time, the entanglement in quantum many- 
body systems, which is the measure of quantum correla¬ 
tions between the two systems, has also become a topic 
of intensive research interest for last several years |1^ 
[ 2 ^. One of the quantities to measure the entanglement 
between the two subsystems is the von Neumann en¬ 
tropy [ 25 I - Q . Consider a bipartite system divided into 
two subsystems A and B of length La and Lb with total 
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length L = La Lb- If the whole system is in a quan- 
turn pure state j'l/') with density matrix p = |'0)(^/)|, then 
the von Neumann entropy S of system A with reduced 
density matrix pA = Avb{p) is defined as 


S = -Tr(pAlogPA)- 


( 2 ) 


The entanglement entropy (EE) or S increases with in¬ 
creasing quantum correlations (entanglement) between 
the two subsystems. EE exhibits distinct scaling rela¬ 
tions at and close to a quantum critical point with the 
shortest length scale of the system. For a critical spin 
chain with periodic boundary conditions where the sub¬ 
system has two boundary points, the entanglement en¬ 
tropy scales as S' = I log La , where c is a universal quan¬ 
tity and given by the central charge of the conformal field 


theory [18|, l20|, 12^ . On the other hand, away from the 


critical point where the correlation length ^ <C La, en¬ 
tanglement entropy is given by S = |log^. The above 
scaling relations are valid for a one-dimensional homo¬ 
geneous system. It is found that some modihcations are 
required in the scaling relations of EE when the system is 
inhomogeneous. Interestingly, in this case also the scal¬ 
ing relations remain same as the homogeneous case with 
a changed prefacto r Cs fF which is called the effective cen¬ 
tral charge [U, [2^-[3J . 

With our understanding of behavior of EE in an equi¬ 
librium system getting better, a considerable amount of 
focus is also given to EE in systems out of equilibrium 
[sB-Ql- The experimental demonstration of such non¬ 
equilibrium dynamics using optical lattices Q also con¬ 
tributes to the sudden upsurge in studies related to deco¬ 
herence and entanglement in out of equilibrium systems. 
One of the ways of generating such a non-equilibrium 
dynamics is a sudden quench. A sudden quench in the 
system can be performed locally or globally. In a global 
quench, a parameter of the Hamiltonian is changed sud¬ 
denly at all the sites resulting to a non-equilibrium dy¬ 
namics. In this process, EE generally shows a linear 
increase in time t up to some time to [13 ■ The local 
quench is defined as a local change of a parameter of 
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the Hamiltonian. For example, the entanglement entropy 
between two critical subsystems A and H of a homoge¬ 
neous one-dimensional chain which are disconnected for 
t < 0 and connected at t = 0 increases as S' = ^ log t for 
t L 113,113; here the final chain is periodic. On the 
other hand, if the final chain is open, the factor 2 in the 
expression of S is not present. Such studies are impor¬ 
tant in the context of information propagation through 
a quantum many body system. 

In this paper, we consider two independent transverse 
field Ising spin chains in the ferromagnetic or critical 
phase. For the times t < 0, the two spin chains which are 
in their respective ground states, are disconnected and 
are later joined at t = 0 (J-quenching). Simultaneously, 
we change the magnitude of the transverse field at a sin¬ 
gle site of one of the chains (/i-quenching). This results 
to a non-equilibrium evolution of the state of the system. 
If we consider only J-quenching of a critical Ising chain, 
both the quantities, namely, LE and EE show periodic 
time evolution [13 ■ However, in this paper, we address 
the effect of /i-quenching along with the J—quenching in 
critical and off-critical systems which results to a non¬ 
trivial evolution of LE and EE when compared to the no 
/i-quenching case. We shall try to understand these re¬ 
sults using the picture of quasiparticles generated due to 
both the local quenches. Our most significant observa¬ 
tion here is the reflection of the quasiparticles at the site 
of /i-quenching. We find some interesting time scales in 
the evolution of both LE and EE when the two quenches 
are performed simultaneously. These timescales can be 
explained successfully using the reflection picture of the 
quasiparticles generated. We have also argued qualita¬ 
tively how the evolution of EE changes as the total sys¬ 
tem moves from deep ferromagnetic phase to critical one. 

The organization of the paper is as follows: we discuss 
the model studied in this paper along with briefly men¬ 
tioning the numerical techniques in section III Al followed 
by a discussion on semiclasscial theory of quasiparticles 
generated in section III HI In section IHII we present the 
results of LE and EE in the critical region for various 
geometries whereas we present the results for the ferro¬ 
magnetic region in section lTVl A comparison between the 
two quenches studied in this paper is made in section El 
using the semiclassical picture. We conclude this paper 
with our main results in section IVTl We have also added 
two appendices at the end of the paper outlining the nu¬ 
merical techniques used in this paper. 


II. MODEL 
A. Exact diagonalization 


where and Jn are the site dependent transverse mag¬ 
netic fields and cooperative interactions, respectively, 
and and are standard Pauli matrices at the lattice 
site n. For the homogeneous case (ft.„ = h and Jn = J), 
the model in Eq. ([3]) has a QCP at J = h separating 
ferromagnetic and quantum paramagnetic phases. Us¬ 
ing Jordan-Wigner transformations followed by Fourier 
transformation for a homogeneous and periodic chain, 
the energy spectrum for the Hamiltonian ([3|) is obtained 
as [isLIiq 


Eg = ±2 J 


\j{h + cosg)2 


J- sin q, 


(4) 


where q is the momentum which takes discrete values 
given by g = 2'nm/L with m = 0 • • ■ L — 1 for a finite 
system of length L. 

On the other hand, such homogeneous systems are very 
rare in nature. One atleast finds some local defects, no 
matter how pure the material is. The general method 
adopted to study systems which are not homogeneous is 
outlined below, which we also use in this paper. Follow¬ 
ing Jordan-Wigner transformation, the Hamiltonian in 
Eq. ([3]) can be described by a quadratic form in terms of 
spinless fermions Ci and c\ 



(5) 


Here, A is a symmetric matrix due to hermicity of H 
and B is an antisymmetric matrix which follows from 
the anticommutation rules of Cj’s. The elements of these 
matrices thus obtained are: 


Ai^j — J- J‘^hi6i^j^ 

^i,j ~ (djJjy-i-i (6) 

The above Hamiltonian can be diagonalized in terms of 
the normal mode spinless Fermi operators rjk given by 
the relation 

Vk = + hk{i)cl), (7) 

where gk{i) and hk{i) are real numbers. In terms of these 
operators the Hamiltonian takes the diagonal form, 

H = J2^k(ilir]k-^\ (8) 

with Afc being the energy of different fermionic modes 
with index k. These A^s are given by the solutions of 
the eigenvalue equations, 

(A-B)(A + B)$fe = Al^k 

(A + B)(A-B)^fc = A^^fc . (9) 


The Hamiltonian we consider here is that of a one¬ 
dimensional Ising chain in a transverse field given by 

H = - + /inO, (3) 

n 


It can be shown that the elements of the eigenvectors 
are related to g and h matrices used to diagonalize 
the Hamiltonian as follows: d>fe(i) = gfe(i) J- hk{i) and 
4'fc(*) = 5fc(*) — hk{i). We calculate LE and EE using $, 
vk, g and h as discussed in appendices lAl andiBl 
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B. Semiclassical theory of quasiparticles 

When a system at zero temperature is taken away 
from its ground state by applying some perturbation, the 
state of the system undergoes a non-equilibrium evolu¬ 
tion with respect to the final Hamiltonian. The initial 
state, which now is an excited state, is a source of quasi¬ 
particles (QPs) corresponding to the final Hamiltonian. 
Recently, such non-equilibrium dynamics have been stud¬ 
ied using a semiclassical picture of quasiparticles gener¬ 
ated for global [13,13 local quenches [3l i where ex¬ 
cellent agreement between the numerics and the semiclas¬ 
sical theory were obtained. We now briefly describe this 
theory of quasiparticles generated which shall be used to 
explain the various timescales observed in our numerical 
calculations. For global quenches from h = 0 to a very 
small h value, it can be shown that these quasiparticles 
are wavepackets of low-lying excitations and discussed 
in details in Ref. [131. Due to the conservation of mo¬ 
mentum, quasiparticles of a given momentum are always 
produced in pairs, the group velocity Vg(k)(= \d£k/dk\) 
of them being equal and opposite to each other. As dis¬ 
cussed in Ref. [13] i these quasiparticles in the small h 
limit can be considered as classical particles (sharply de¬ 
fined QPs) which when crosses a site, simply flips the 
spin at that site. Though this picture is discussed for a 
very specific quench (a small quench), it has been veri¬ 
fied for stronger quenches and also for quenches in the 
paramagnetic phase with slight modifications. It is also 
argued that these quasiparticles are no longer point par¬ 
ticles, but are extended objects as the critical point is 
approached due to large correlation length. In the fol¬ 
lowing sections, we shall try to explain our numerical 
results atleast qualitatively with this point like picture 
of quasiparticles for spin chains in the ferromagnetic as 
well as in the critical region. 


III. LOSCHMIDT ECHO AND 
ENTANGLEMENT ENTROPY FOR CRITICAL 
CHAIN 

We first study double quenches for a critical chain 
where already some work has been done in Ref. [ 13 , 113 | 
for local J-quenches. As discussed before, we consider 
simultaneous application of two types of local pertur¬ 
bations to the system and study the time evolution 
of Loschmidt echo (LE) and von Neumann entangle¬ 
ment entropy (EE) as a result of such quenches. Ini¬ 
tially the spin chain is prepared in the ground state of 
H = Hi + H 2 where Hi and H 2 are the Hamiltoni¬ 
ans of two decoupled homogeneous Ising chains of length 
Li and L 2 , respectively, with open boundary conditions 
[Jli = Jl 2 = 0)- Two simultaneous quenches are per¬ 
formed at t = 0, namely, (i) the two spin-1/2 chains are 
suddenly connected together resulting to a chain of total 
length L = Li + L 2 , and, (ii) the transverse field at a 
particular site L' belonging to either the chain 1 or 2 is 


changed from h to h + 6. The system then evolves with 
the final Hamiltonian 

Hf = Hi+H 2 + Hi^ - Sal ,, ( 10 ) 

where HI2 defines the connection between the two spin 
chains of length Li and L2 and is of the form 
At the same time, the term in the Hamiltonian 

corresponds to the ft,—quenching which changes the mag¬ 
nitude of transverse field at site H from h to h + 6. We 
incorporate these quenches numerically by considering a 
single spin chain of total length L = (Li + L 2 ) with the 
first Li spins forming the system 1 and the remaining 
system 2. They are disconnected at t < 0 by putting 
Jli = 0 which at t = 0 is then increased to J, also the 
interaction strength at all the other sites. Eor all our 
calculations, we have set J = 1. The details of the nu¬ 
merical calculations for LE and EE are outlined in the 
appendix, see also Refs. [IS [ 13 . 

As we switch on the two local perturbations discussed 
above, there is a local increase in e nerg y of the system at 
the site of local perturbations [^ [3l ■ These sites then 
become the source of quasiparticle production. Hence¬ 
forth, we shall call the quasiparticles created due to the 
ft-quenching at L' as QP^ , and the corresponding left 
and right moving quasiparticles as QPl and QP}^, respec¬ 
tively. Similarly, the left and right moving quasipartciles 
created at the site Li of J—quenching shall be called as 
QP/ and QP^- Below, we present our results for the evo¬ 
lution of LE and EE for various cases or geometry and 
discuss these results in the light of quasiparticles propa¬ 
gating in the system. 



FIG. 1: The plot shows LE as a function of time for different 
values of <5 with J—quench at Li — Z//2 and the ft—quench 
at site L' = L/S. The transverse field at L' is changed from 
1 to 1 + 5 . For the J quenching alone (i.e., <5 = 0 case), the 
LE shows peak at ts = LjVina.x = 150 and T = 21^ where 
Umax = 2 and L = 300. By applying two local perturbations 
simultaneously at time t = 0, we observe a small peak at t' = 
50 and comparatively a stronger peak at t” = 200. We also 
note small fluctuations near ti = 100 which is more clearly 
seen for <5 = 1 curve. 
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A. LE and EE for Li = L2 


In this section, we consider J—quenching at Li = L 2 = 
L/2 and the h-quenching at some site L' of the total 
spin chain of length L at time t = 0. Let us first dis¬ 
cuss the dynamics of LE. In general, for the better read¬ 
ability of the paper, we shall assume L' < Li < La 
(explained later in the context of EE) through out the 
paper but the case with L' > Li is also presented and 
discussed in the caption of various figures. Due to the 
propagation of the generated quasiparticles, we expect 
four time scales which are the times of come back of 
the quasi-particles at the source point after getting re¬ 
flected from boundaries of the chain, thus removing the 
effect of their dynamics. When this happens, the over¬ 
lap in the definition of LE increases and shows a peak. 
These time scales are determined by the fastest moving 
QPs with maximal group velocity r^max = max^ Vg{k) 
and are given by ti = (2L')/vmax (time of come back 
of QPl), t 2 = 2(L — L')/vms,x (time of come back of 
QPk)^ ^3 = (2Li)/uniax (time of come back of QPl) and 
t 4 = (2L2)/fmax (time of come back of QP^)- With 
Li = L 2 , the values of and ti are equal, but it is not the 
case in general. Other than the above mentioned obvious 
time scales, two more time scales are observed numeri¬ 
cally, t' given by 2|Li —L'l/umax, which appears to be the 
time taken by QP]^ (for Li > L') to reach L' and get par¬ 
tially reflected at L' where it finds a change in potential 
from h to h-\-5. The second time scale is same as t 2 , but it 
is due to QP^ and is given by t" = 2{L —L')/vmax, which 
is the time taken by the {QPr) to get reflected at 
L' (right boundary) and come back to Li after getting 
fully (partially) reflected at the right boundary {L'). It 
is to be noted that Umax for the homogeneous transverse 
Ising model with elements as defined in Eq. [5] is 2 at the 
critical point and also in the paramagnetic phase. We 
shall use the same value of Umax in our case also since the 
numerically obtained value of Vmax by differentiating the 
eigenvalues is also close to 2. 

We show the evolution of LE in Fig.[T]for different val¬ 
ues of 6. Let us first concentrate on the results of LE 
with i5 = 0. Since LE is the overlap of two wavefunc- 
tions which had unity overlap initially, it starts decreas¬ 
ing from one at t = 0 till both the quasiparticles reach 
the boundary at t = LjA. (since L\ = L 2 = L12 and 
VmaK = 2), where it gets reflected. Intuitively, during 
its return path, QPs will undo their effect of dynamics. 
Thus, we expect to see a decrease in LE till the reflec¬ 
tion of the first quasiparticle {i.e., till t=L/4) after which 
there is an increase till the quasiparticle reaches its origin 
or till t = L/2{= 2Li/2). After this time, the initially 
left (right) moving quasiparticle will move to the sys¬ 
tem 2 (system 1) and eventually come back to its origin 
showing a peak at t = T = L, which is also the time pe¬ 
riod of the quasiparticles, see Fig. [T] When the second 
quench or the h-quench is also performed simultaneously, 
we expect to see some structures close to the time scales 
mentioned before, i.e., at ti,t 2 ,t', t" along with the (5 = 0 



FIG. 2: The plot shows LE as a function of time when h- 
quenching of strength 5 = 1.0 is performed at different sites L' 
of the total chain with J-quenching fixed at Li = L/2. Here 
L' = 0 corresponds to the case of J-quenching alone where 
ts = ti = 150. For L' = L/3 ,t' = 50 and t" — 200 whereas 
t' = 90 and t" — 240 for L' = L/5. All these timescales 
are clearly seen in the above figure. We also observe small 
perturbations at which is not very clear. 



FIG. 3: Entanglement entropy as a function of time after a lo¬ 
cal J-quench in the middle of the chain (Li = L/2) along with 
/i-quenching at L' = Lj/t for different interaction strengths 5. 
Here, we consider La = Li = L 12 as the subsystem with total 
system size L = 300. The time scales t'/2, t' and t” can be 
seen in this figure which agrees well with our explanations. In 
this case, t' = 75 and t" = 225. Another time scale observed 
is around t = 187 which is when QPji enters system A result¬ 
ing to an increase in EE as QPl is in system B during that 
time. Note that the effect of QP^ is very small compared to 
QP^. 


structure. The numerical results of such double quenches 
are shown in Fig. [T]for various coupling strengths 8 and 
fixed L'. As expected, the decay in LE is stronger with 
increasing strength of the local h-quench or 8. On the 
other hand. Fig. [2] shows the variation of LE(t) when the 
/i-quench is performed at different positions of the chain 
with fixed (5=1 demonstrating the above mentioned time 
scales more clearly, especially the variation of t" with L'. 
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FIG. 4: Time evolution of entanglement entropy after a local 
J-quench with h-quenching at different sites L' for 5 = 1.0. 
For L' = L/& = 50, the deviation from the no h-quench case 
starts at t' /2 = 50. A sharp decrease is observed at t = 
100 when QPi, gets reflected at L' and returns to Li = La- 
There is a sudden increase in EE after t = 150 when QP]i 
enters system A whereas QPl is still in system B. We again 
see a dip at t" = 250 when QPl and QPr exchange their 
systems. Similarly, for the case L' — 2L/3 = 200, QP^ enters 
system A at t = t' /2 = 25 when the deviation from the single 
quench case appears. We see a sudden decrease of EE a.tt = t' 
when QPr enters system A after getting reflected at L' so that 
both the quasiparticles are in system A. After t = 150, once 
again EE increases as QP^ enters system B. A sharp dip is 
seen at t = t” = 200 when QP^ and QP]i exchange systems. 

We find good agreement between the timescales proposed 
above and the numerics, see the caption for more details. 

In conclusion, we find that the dominant peaks are due 
to QP^ at times t', t", and ^ 4 , where ts = ^4 in this case. 
We also note that the peak at t" is a strong peak which 
may be due to the fact that at this time three different 
quasiparticles return to their origin after reflections at 
various points as discussed below: (i) QP}i after reflection 
from the right boundary, {\\)QP^ after reflection from 
right boundary and a second reflection at L' causing it to 
return to Li (iii) QPf^ after reflection at L' and a second 
reflection at right boundary resulting to its return to Li. 
Finally, a peak is also observed a.t t = L which is the 
return time of all the fastest quasiparticles back to their 
origin when there is no reflection at L', thus giving us a 
hint that there may be a transmitted component of the 
quasiparticle also. We shall comment more on it after 
discussing the results in the ferromagnetic phase. 

It is to be noted that the presence of time scale ti due 
to QP^ is almost negligible in these figures, though we 
do observe some perturbation at this time. On the other 
hand, it is too early to discard the presence of QP^ as 
its effect is very clearly observed in the evolution of EE 
as explained in the next paragraph, thus ruling out the 
possibility of absence of such quasiparticles. We shall try 
to argue about the absence of scale in the evolution of 
LE in section 0 

We now focus on the entanglement entropy as a func¬ 


tion of time for the above scenario. In this case, another 
parameter is the size La of system A of which we cal¬ 
culate the entanglement entropy with the remaining sys¬ 
tem of size L — La- Let us first consider the simplest 
case where Li = La, i.e., the location of J-quenching 
also determines the size La of system A. Interestingly, 
the bipartite EE of two critical transverse Ising chains 
can also detect the response of h-quenching (see Eig. |31 
a. We see that for J quenching alone or for the sin¬ 
gle quench, the entanglement entropy shows perfect pe¬ 
riodic oscillations with dips at ^ 3=^4 = 2 Li/umax, also 
discussed in Ref.ji^ using conformal field theory. This 
can also be explained using the quasiparticle picture. A 
pair of quasiparticle will increase the entanglement be¬ 
tween the system A and the rest if one of them is in 
system A and the other is in B. The quasiparticle pairs 
are generated at Li = La = Lj^ and travel in oppo¬ 
site directions resulting to an immediate increase in S'(t) 
for t > 0. This is not the case when Li ^ La and will 
soon be discussed separately. As both of them reaches 
the boundary at t = Lj^ and gets reflected, S starts 
decreasing and eventually shows a dip when t = L/2 
after which both the quasiparticles belonging to a pair 
exchange their systems and once again EE increases af¬ 
ter t = L/2. At T = 2L/umax(= L), both the quasi¬ 
particles arrive at the starting point and S shows a dip 
once again after which the pattern repeats. This is also 
shown in Eig |4] with L' = 0.0. If we now perform local 
h-quenching at a general site L' of the spin chain, we ob¬ 
serve that the evolution of EE in double quenches follows 
the single quench case but accompanied by deviations at 
certain times which can once again be explained using 
the quasiparticle picture. We observe that the double 
quench case follows the single quench case (<5 = 0.0) till 
t = t '/2 after which there is a sudden deviation or in¬ 
crease from the single quench case. This is because one 
of the quasiparticles {QP}^) produced at L' (which is not 
present in single quench case) enters the system B at this 
time whereas the other quasiparticle of the same pair re¬ 
mains in system A. This results to an extra increase in 
S- On the other hand, at t = t', the QPj^ reaches back 
to Li = La after reflection at L' where we find a sharp 
decrease in S as both QPf^ and QP^ are now in system 
B. The EE keeps on decreasing (after a slight increase) 
till QPji enters system A at t = 2(L — Li)/umax as a 
result of getting reflected from the right boundary, after 
which we observe a sharp increase in EE. One more time 
scale observed corresponds to t" = 2(L — L')/umax which 
was also observed in LE. At this time, the QP]i and QPl 
return back to La and exchange their systems, as also 
discussed before with reference to LE. 

The main difference between the analysis of LE and 
EE is that in the case of LE, we are interested in time 
scales at which the produced quasiparticles come back 
to its origin. In the case of EE, we are interested in the 
timescales in which one of the quasiparticles belonging to 
a pair crosses one system and goes to the other system. 
If this crossover results to both the QPs to be in the 
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same system, then EE decreases, otherwise it increases 
Here, since we considered the geometry where Li = L 2 = 
La, many of the time scales are hidden. In the next 
section, we apply the same ideas to the case when Li ^ 
L 2 but Li = La along with a special discussion for the 
most general case when Li ^ L 2 ^ La, and verify the 
quasiparticle picture proposed. 



FIG. 5: The plot shows LE as a function of time for double 
quenches with L\ — 100, La = 100 and L = 300 and different 
L'. The first peak of LE occurs at time t' — 2(Li — L')/vniax 
(t' = 25 and 50 for L' — L/A and L/6 respectively). The 
other time scales are fa = 100, ti — 200 and t". We note that 
t” = 225 and 250 for L' = 75 and 50, respectively. 


B. LE and EE for Li yf L2 

We are now interested in the local quenching of asym¬ 
metric spin chain (Li fy ^ 2 )- Let us first study the evolu¬ 
tion of Loschmidt echo. Eor J-quenching alone [1^, the 
LE shows three time scales given by fy, fy and the time 
period T discussed in section IIII Al (see L' = 0 plot of 
Fig.[5|). The nature of these plots are discussed in details 
in Ref. [i^ using conformal field theory. When the trans¬ 
verse field term at L' is changed from h to h + 5 together 
with J-quenching, we expect to see the following addi¬ 
tional time scales in parallel with the discussion in the 
previous section: ti = (2L')/?;max, t' = 2{Li - L')/^max 
, t" = 2{L — L')/'Ui„ax Fig. [5] shows LE as a function of 
time after J-quenching at Li = Lj^t and h-quenching at 
different sites L' . All the above mentioned time scales 
can be clearly seen in this figure except fy. We shall try 
to argue for this latter. 

Let us now move on to the calculation of EE for the 
same situation {Li ^ L2, but Li = La)- Fig.|6]shows the 
time evolution of EE after single and double quenching 
at time t = 0. One can observe the difference in LE and 
EE between the times ts and t 4 for J = 0 case (see Fig. [5] 
and Fig. [6]). In this time range LE remains constant. 
On the other hand EE decreases slowly and then starts 
increasing at t = t 4 as discussed in Ref. [l^ , the increase 
being due to arrival of QP^ in system A. Moving to the 


double quenches, the discussion is almost same as for the 
case Li = L 2 = L/2 in the previous section. EE more or 
less follows the J = 0 case and we see special time scales 
at t'/2, t' and t". See caption of Fig. [5] for more details. 



FIG. 6: Time evolution of entanglement entropy for the same 
case as in Fig. [S]but different L'. For L' = L/4, the deviation 
from single J-quench case appears at t' 12 = 12.5 whereas at 
t' = 25, QPI enters system B after getting reflected at L' 
where its other partner QP^ is already present. The decrease 
in EE continues till QP^ enters system A at t — L2 = 200. We 
also see a dip at t = t” = 225 where QP^ and QPr exchanges 
their systems. Similarly, one can argue for the evolution of 
EE when L' = -tLjA:. The deviation from the single quench 
case begins at t — 62.5. In this case, QPl enters system B 
at t — 100 which causes a sharp decrease at this time. On 
the other hand, QP^ enters system A at t = 125 resulting 
to an increase in EE as its other counterpart is still in B. 
The natural increase at t — 200 which is there for only J- 
quenching case can also be observed. This might be due to 
the fact that any reflection at L' is not perfect and there is 
a possibility of getting a transmitted component of the QP 
wave, also discussed in sections [TV] and ED The time scales 
t” — 225 and T — 300 are also present. 

We now consider the most general situation with Li fy 
L 2 fy La and calculate the time evolution of EE. Here, 
J-quenching is again performed at Li, but the subsys¬ 
tem is assumed to be of length La = L/2, different 
from Li- The EE as a function of time after the double 
quenches is shown in Fig. [T] Let us define l{= La — Li) 
as the distance between the right end of the subsystem 
A and the site of J-quenching. As in other cases, we 
discuss explicitly the case with L' < L4 < La below, 
and try to present some other examples through the fig¬ 
ures. For J = 0, EE remains at a constant value for 
small times and starts increasing at t = Z/umax when 
QP^ hits at the boundary of the two subsystems and 
enters the subsystem B [s^. Note the contrast be¬ 
tween Li = La where S(t) increases immediately and 
Li fy La where S{t) is constant initially. The EE shows 
a sharp decrease at t = (2Li -|- l)/vmax when QP^ en¬ 
ters system B after getting reflected from the left bound¬ 
ary. Similar to the previous case (i.e., for Li fy L 2 and 
Li = La), there is a ’decay region’ between the time 
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FIG. 7: Time evolution of entanglement entropy after a local 
J-quenching at L\ = L/S and /i-quenctiing at different sites 
L' with (5 = 1.0. Here, the subsystem is of length La ~ 
I//2 which does not coincide with cut resulting to few more 
relevant time scales. For L' — L/5 = 60, the first deviation 
(or increase) from the single quench case appears at t = 45. 
The EE decreases at t = t' = 65 when QP^ and QPr are 
in the same subsystem B. The next increase in EE would be 
at t = 175 when QPr enters subsystem A. The split of time 
scale t" , as discussed in the text, can also be seen with dips 
at i!{ = 215 and = 265. The case with L' = 3Z//4 is all the 
more interesting. The deviation occurs a.t t = {L' — La)12 = 
37.5. The sharp decrease in this case occurs at t = 100 which 
is the time taken by QPr to get reflected at L' and enter 
system A so that both QPr and QPr are in system A. But 
at t = 125, QPr enters system B after getting reflected from 
the left boundary resulting to an increase in EE. t" = 200 
(due to QPr) and t 2 = 250 (due to QPr). We note extra 
dips around t = 187 which seems to be due to QPr entering 
system B after reflection from the left boundary. 


range [(2Li + 0/t^max, {L 2 + 7.B)/t^max] where EE decays 
very slowly, after which there is an increase as both the 
QPs are now in different subsystems. In this case, since 
the site for J-quenching does not coincide with the sub¬ 
system size, there are many additional time scales as dis¬ 
cussed below, the appearance of which puts the picture of 
travelling quasiparticle on stronger footing.. Let us now 
come back to the double quenches. Following the double 
quenches, EE more or less follows the single quench case. 
The first deviation resulting to an increase in EE (similar 
to the one discussed before) occurs at t = (L^i —L')/uniax, 
when QP^ enters the subsystem B. By definition, the 
time scale t' is the time taken by the QPr to get re¬ 
flected at L' and come back to system A, which in this 
particular case is given by t' = (2(Li — L') -|- O/i'max- 
On the other hand, the time scale t", defined as the time 
taken by QPr (or QP^ to undergo double reflections 
at L' and one of the boundaries, gets divided into two 
scales. This is because the distance travelled by QPr is 
smaller than QP^, which was not the case in our previ¬ 
ous discussions where Li = La- Let us define these two 
timescales by t'/ = {2{Li—L') + 1 + 2Lb)/‘2 (for QPl) and 
<2 = t" -I- 2l/vms,x (for QPb)- All these time scales are 


clearly shown in Fig. [71 We would like to point out here 
that the basic physics related to tracking of QPs remain 
same when we change the position of L' , but the formula 
for these time scales may have to be changed as can be 
seen in the L' = 3L/4 case discussed in Fig. jT] Also, the 
above discussion will be correct if [La — L')/2 < Li— L', 
i.e., QPb reaches La before QPr. In the opposite case 
also, one needs to simply apply the same ideas to get the 
right picture of dynamics. It is also to be mentioned that 
we do observe some extra time scales, some of which can 
be explained and are discussed in the caption of Fig. jT] 


IV. ENTANGLEMENT ENTROPY FOR A 
FERROMGANETIC CHAIN 

In this section, we briefly discuss the evolution of en¬ 
tanglement entropy when the total spin chain is in the 
ferromagnetic phase. Here, we consider local quenching 
of asymmetric spin chains (Li 7 ^ L 2 ) with Li = La. We 
concentrate upon two different cases to calculate EE af¬ 
ter single or double quenches : one where the total spin 
chain is deep inside the ferromagnetic phase (see Fig. [ 5 ]) 
and the other where the spin chain is close to the critical 
point (see Fig. HD. Let us first consider the spin chain 
with h = 0.5 at all sites. Fig. |5| shows time evolution of 
EE after single quench at Li {J quench) and also after the 
double quenches, namely, h quench at L' and J quench 
at Li. For the single quench {U = 0), EE detects Q and 
t 4 successfully. Note the difference between the critical 
and ferromagnetic region for times up to Q. In the criti¬ 
cal case, the EE increases at t = 0 followed by a decrease 
which starts around Q/2 when the QPr gets reflected 
from the boundary, though the decrease is sharper at Q. 
On the other hand, in the ferromagnetic region, we see 
a sudden increase in EE followed by an almost constant 
EE region up to Q (no decrease at Q/2) after which it 
decreases suddenly. This hints to the fact that in the 
ferromagnetic region, QPs are more point like particles, 
and hence its location can be known precisely. But in 
the critical case, these QPs are extended wavepackets as 
also mentioned in Ref. [^, and hence the reflection at 
the boundary is felt also at Li. Similar to the critical 
case (see Sec. IIIIBI) . EE decays between times Q to t^. 
In this time range the fastest moving quasi-particles do 
not contribute in the EE. Let us now discuss the time 
evolution of EE after double quenches. One can observe 
clearly the time scales t'/2 and t' from Fig. [51 Interest¬ 
ingly, EE starts decreasing after t' and it continues up to 
t 4 . The sharp increase in EE at is due to the fact that 
at O, QPb enters system A whereas QPr is still in sys¬ 
tem B. It is to be mentioned that in the ferromagnetic 
case, t" is not clearly visible, which once again can be 
explained due to the point like nature of QPs in the fer¬ 
romagnetic region. For t < t", QPr is in system B and 
QPb is in system A. They exchange their systems at t", 
thus contributing to the entropy equally for t < t" and 
t > t". On the other hand, the critical case distinguishes 











between QPs approaching Li and moving away from Li 
due to the finite extent of QP. 

We now move to explain the quenching results when 
the final system remains close to the quantum critical 
point. The evolution of EE as a function of time fol¬ 
lowing single and double quenches is shown in Fig. [5] 
for h = 0.99. In this case also, we observe a sudden 
increase/deviation from the single quench case at t' 12 
followed by a sharp decrease at t'. We also observe an 
increase immediately after t' which is different from the 
h = 0.5 case and similar to the critical case. This may 
be because the QP which is now more like an extended 
object with extended wavefunction is only partially re¬ 
flected at L' as compared to localized QP deep inside the 
ferromagnetic phase having less wavelike properties. Rest 
of the discussion is the same in this case and discussed 
in details in the caption of Fig. [9] 



FIG. 8: The plot shows EE as a function of time for single 
and double quenches when the whole chain is non-critical {h — 
0.5) with Li = 100, La = 100 and L = 300 and different L'. 
For L' = i/4, the deviation in EE from single quenching 
case starts at t = t'/2 where t' = 50, with Umax = 2h = 1. 
For t > 50, both QP^ and QP^ are in system B leading to 
decrease in EE which continues up to ti = 400. Similarly, 
for L' = L/5 one can find t' = 80 and the figure shows the 
expected behavior. The absence of t” is explained in the text. 


V. COMPARISON BETWEEN QP^ AND gP^ 

In this section we try to argue why the effect of QP^ 
is stronger than QP^. As discussed before, the initial 
state, which is no longer the ground state of the final 
Hamiltonian, is a source of quasiparticles. Also, since 
the perturbation studied in this paper is local and very 
small, the quasiparticles are produced at the site of per¬ 
turbation only, i.e., at L' and Li. These quasiparticles 
have energies given by the eigenvalues of the final Hamil¬ 
tonian. Let k identifies the quasiparticle rjk having eigen 
energy Afc. Quasiparticles of energy A^ are produced at 
the perturbation site with probability fk- Numerically, 



FIG. 9: EE as a function of time for the same situation as in 
Fi^but fixing h at 0.99. In this case the value of Umax is 1.98. 
This gives t' = 25.25 for L' = L/4. We do see timescales t'/2 
and t' along with ts 101) and (~ 202). We also observe 
a peak near t" (~ 227). 

one can obtain fk by calculating the expectation value 

fk = {fpi\vlvk\tpi) 

which is proportional to the number of quasiparticles r]k 
present in the initial state |'0i)- This expression can be 
written in terms of Tfc of the final Hamiltonian and 
the matrix G® (see Appendix IbI) with respect to the ini¬ 
tial Hamiltonian. A comparison of fk for the h-quench 
alone {h changed from 1 to 2 at L'), J-quench alone (0 
to 1 at Li) and both quenches together is shown in Fig. 
1101 Clearly, quasiparticle creation probability is an or¬ 
der of magnitude higher in case of J-quench alone when 
compared to h-quench. This hints to the fact that the 
J-quench is the main source of quasiparticle production 
and hence our numerical results are dominated by the 
dynamics of QP^. 



FIG. 10: Variation of fk with k for a critical chain. The 
dashed line or the lowest curve corresponds to the h-quench 
case which clearly is an order of magnitude smaller than the 
J-quench alone. Similar behavior is also observed for a ferro¬ 
magnetic chain. 
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VI. CONCLUSIONS AND DISCUSSIONS 


In this paper, we studied the effect of two simultaneous 
local quenches in an otherwise uniform transverse Ising 
chain of length L with open boundary conditions. Ini¬ 
tially, the system is prepared in the ground state of the 
transverse Ising chain having a uniform transverse field h 
and the interaction strength set to unity at all sites except 
Jli and Jl where it is zero. The first quench corresponds 
to sudden increase of Jl^ from zero to I, and the second 
quench involves the sudden change of the transverse held 
from hio h + 5 aX site L'. We argued that the sites of the 
two local quenches are source of quasiparticle production 
as there is a local increase in energy due to the quenches. 
These QPs are wavepackets of low lying excitations of 
the hnal Hamiltonian. As discussed in Refs [13,13 j the 
QPs are localized in the ferromagnetic region and behave 
more like classical particles, whereas they are extended 
objects/wavepackets as the critical point is approached. 

We numerically studied the evolution of Loschmidt 
echo and entanglement entropy after the double quenches 
and explained the evolution using the quasiparticle pic¬ 
ture. The envelope of the curve is dictated by the fastest 
moving quasiparticles. We showed taking examples that 
most of the timescales can be explained using the prop¬ 
agation of quasiparticles. The most interesting phenom¬ 
ena that is observed numerically is the partial or full re¬ 
flection of the quasiparticles at L', the site of ft-—quench. 
Only if we include such a phenomena that we can explain 
certain numerically observed time scales. As mentioned 
in details in the paper, the most relevant time scales are 
t', ta, ti and t" (for Li = La) which are all due to the 
QP^ pair, or the pair produced at the J-quenching site. 
The presence of the other set of quasiparticles, namely, 
QP^ is clearly seen in the evolution of EE. We have shown 
that the probability of quasiparticles produced due to h- 
quench is roughly an order of magnitude smaller than the 
J-quench which could be the reason for stronger effect of 
QP^ in the evolution of LE and EE. 

The double quenches deep inside the ferromagnetic 
phase can be very nicely described by the point like quasi¬ 
particles where all the time scales are sharply observed. 
We have contrasted this ferromagnetic case with the dou¬ 
ble quenches in the critical phase and proposed the rea¬ 
sons for their differences. It seems that the reflection 
of QP^ at L' in a critical chain is only partial having a 
transmitted component also. This can be attributed to 
extended wavepacket nature of quasiparticles at the crit¬ 
ical point. One can then explain the decrease of EE at 
t = ts, slight increase of EE for t > t' after the sharp 
decrease at t', and the dip at t = T. 

Our main aim in this paper is to study the dynamical 
evolution of LE and EE after double quenches and see 
if one can explain the behavior, atleast qualitatively, us¬ 
ing propagation of quasiparticles. We have demonstrated 
here that this indeed is possible. Though we can not pro¬ 
pose a general formula for all the timescales involved as 
it depends on which quasiparticle arrives at the subsys¬ 


tem first, which in turn depends on the location of L', 
Li and La, but the basic idea gives us the right picture. 
We have checked this for other cases also which are not 
presented in this paper. The quasiparticle picture does 
explain many features, if not all, of the dynamical evolu¬ 
tion of LE and EE that occur in double quenches studied 
here. We have provided some arguments for the imme¬ 
diate increase of EE after t' for a critical chain which is 
related to the extended nature of the QPs at the critical 
point. 
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Appendix A: Loschmidt echo for a general quadratic 
fermionic system 


Here we shall discuss the method for evaluating the 
time evolution of LE in real space for a general quadratic 
fermionic system 0- We rewrite the Hamiltonian in 
Eq. ([5]) in the following form 

H = ic'fHC, (Al) 

where = ^c|,... c^, ci,... cj^ and LL = ® A + 

0 B. 

The Loschmidt echo, defined in Eq. o, can be evalu¬ 
ated in this case by the following relation 00 

£(t) = |(V^(0)|e-*‘"QV'(0))| = |det(l -R + Re--«^*)|, 

(A2) 

where Hf is the final Hamiltonian after double quenches. 
R is the 2L x 2L correlation matrix whose elements 
are two-point correlation functions of fermionic operators 
Rij = {^|;{0)\C}Cj\^|:{0)), where |V’(0)) is the ground state 
of the initial Hamiltonian Lli. Following some mathe¬ 
matical steps we can the matrix R in terms of g and h 
matrices (see discussion around Eq. [7|) 


R = 


h' (h*)T 
gi (h‘)T 


h* (g')'^ \ 

g‘ (g')’^ J ' 


(A3) 


Here, the g* and h* matrices are calculated relative to the 
initial Hamiltonian Hi. For more details, see Ref. 0 


Appendix B: Time dependant entanglement entropy 
for a general quadratic fermionic system 

We outline here steps for calculating the time evolution 
of EE after sudden quenches [H, [5l[- As discussed in 
the text, our model is reduced in quadratic form using 
spinless fermions. Let us define two Clifford operators 
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which are also related to the Majorana fermion operators 
a 2 i-i and a 2 i in the following way, 

Ai = c} + Ci = a 2 i-i and Bi = c] - Ci = ia 2 i. (Bl) 

The operators A and B can be written in terms of the 
free-fermion operators r]k 


by an orthogonal matrix say, V. Therefore the eigenval¬ 
ues of T"^ are purely imaginary of the form it'; with 
I = 1,2,-■■ ,La- This can be used to write the re¬ 
duced density matrix as a direct product of La uncorre¬ 
lated modes pA = where each gi has eigenvalues 

(1 ± vi)l2. Thus the bipartite EE for pA is the sum of 
entropies of La uncorrelated modes given by 


L L 

A = Bi = '^'iik{i){vl-Vk)- (B2) 

k —1 k —1 

The time evolution of these operators are obtained 
from the time dependence of fermionic operators, i.e., 
r]k{t) = where pk and Ak are quasiparticles and 

eigen energies corresponding to the final Hamiltonian Hf. 
This is given by 


L 


Ai{t) — [{AiAj)tAj -l- {AiBj)tBj \, 
1=1 
L 

BM = E A + , 

1=1 


(B3) 


Sl{La) 


;=i ' 


l-Vl 

2 



The time-dependent expectation values in Eq. 
calculated using Eq. (IB3I) as 


(B6) 

are 


{a2l{t)a2m{t)) = -{Bl{t)Bm{t)), 
{a2l-l{t)a2m-li.t)) = -{Al{t)Amit)), 
{a2lit)a2m-l{t)) = -i{Bl{t)Am{t)), 
{a2i-i{t)a2m{t)) = -i{Ai{t)Bm{t)) ■ (B7) 


Thus we shall get the time-dependent correlation matrix 
T"^ and for each time we can obtain EE from Eq. (IB6I) . 
The elements of the matrix are given by 


where the time-dependent contractions are 

L 

{A^Aj)t = y^COs(Afct)^/e(z)^Ai(j) , 

k^l 

L 

{A^BJ)t = {BjAi)t=i^sm{Akt)^k{i)^kij), 

k^l 

L 

= ^cos{Akt)'^k{i)'^k{j) - (B4) 

k^l 

The total system is divided into two subsystems A and 
B of length La and Lb respectively. To evaluate EE 
between these two subsystems after the double quench, 
we have to calculate the reduced density matrix pA of 
the subsystem A. This can be reconstructed from the 
2La X ‘2La correlation matrix of the Majorana operators 


jam(^)an (t) — ^m,n A ^(E )mni (B5) 


B2i-l,2m-l “ ^ E! ^kTk,{'^lBki)t{AmAk2)t 

ki ,k2 

+ * E Gl^kAAiAkMAmBk2)t 

ki ,k2 

^ 21 - 1 , 2 m = Y. (^k2kAAlAkMBmBk2)t 
ki ,k2 

- Y (^Ak2iAiBkAABmAk2)t 

ki ,k2 

= -Y^^U^iBlBkAAAmAk^h 

kik2 

+ Y (^U^iBiAkAAAmBk^h 

ki ,k2 

^hm = -^Y B^AkABlAkAt{BmBk2)t 

ki ,k2 

+ lYGAk^iBiBkAABmAk^h- (B8) 

fcl ,fe2 


where m, n = 1,2, S,--- , 2 La and lipi) is the initial 
ground state. The matrix E'^ is an antisymmetric ma¬ 
trix which can be brought into the block-diagonal form 


where = —Yhk^ki^A^kiAA is the equilibrium 

correlation Function which is calculated with the initial 
Hamiltonian Hi. 
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